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Abstract 

The scaling of nuclear spin-dependent parity violating effects with increasing nuclear charge Z 
is discussed in two series of isovalent open-shell diatomic molecules. The parameter W a character- 
ising the strength of parity violation in diatomic molecules is calculated in the framework of the 
zeroth-order regular approximation (ZORA) and found to be in good agreement with the R{Z)Z k 
scaling law derived for atoms in which R(Z) represents a relativistic enhancement factor. The in- 
fluence of electron correlation is studied on the molecular level, with spin-polarisation effects being 
conveniently accounted for by a previously established approximate relation between the hyperfine 
coupling tensor and W a . For high accuracy predictions of parity violating effects in radium fluoride 
the necessity for systematically improvable correlation calculations is emphasised. 
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INTRODUCTION 



Molecular properties depending strongly on the behaviour of the electronic wavefunction 
in the vicinity of the nuclei display a pronounced dependence on the nuclear charge Z . This 
is known for a long time and often employed in atomic physics to obtain qualitative estimates 
(see e.g. Ref. [1]). For systems with more complicated electronic structure, in particular for 
molecules bearing nuclei with various Z values, a simple Z-scaling of such molecular prop- 
erties is not a priori guaranteed, as molecular properties are also depending on the specific 
nuclear arrangement, which could in principle hamper a direct comparison (see e.g. discus- 
sion in Ref. 2J). Establishing general scaling laws also for complex molecules would present 
a great advantage, as scaling laws can be used for inexpensive order-of-magnitude estimates. 
Quantum chemical calculations allow to scrutinise proposed scaling laws for a given property 
and investigate the form of the dependence on Z . This has been applied in some detail in 
Ref. [3[ for nuclear spin-independent parity violating effects in chiral molecules containing 
atoms from various rows of the periodic table (see Refs. 4|-|6) for reviews on molecular parity 
violation). In the present article we are calculating nuclear spin-dependent parity violation 
interaction in alkaline earth metal monofluorides (Mg-Ra)F (group II monofluorides) and 
(Zn-Cn)H (group XII monohydrides) as an example for open-shell systems. Results and 
conclusions presented below were reported by the authors on several workshops and confer- 
ences during the years 2010 and 2011 and explicitly foreshadowed in Ref. [7J. A recent paper 
took up the idea that we reported on those occasions, which motivates us to present here 
our results and to comment, in particular, on the inclusion of electron correlation effects in 
molecular systems. 



NUCLEAR SPIN-DEPENDENT PARITY VIOLATION 

One of the properties that is predicted to depend heavily on Z is the nuclear spin- 
dependent parity-odd (P-odd ) interaction [MST> — W ), whose effective operator looks in 

a four- component (relativistic) framework like [8j: 

G x -. 

Kv = -m Yl kA > A d ' Ia paW, w 

* A,i 

where Gp is Fermi's constant of the weak interaction, k^A is an effective parameter describ- 
ing MST> — VV interactions for nucleus A (caused both by the nuclear anapole moment and 
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by weak electron-nucleon interactions, see Ref. 
density distribution of nucleus A, respectively. 



[91]), 7^4 and pa are the spin and nuclear spin 
For a. one uses 



a 




(2) 



with <r being a vector of the 2x2 Pauli spin matrices o~ x , o ~ y , <j 7 . and with being a 2 x 2 zero 



matrix. The anapole moment was proposed by Zeldovich lCj, quickly after the discovery 



of parity violation in processes mediated by the fundamental weak interaction. The nuclear 
anapole moment [3] has received great interest in atomic and molecular physics as it is 
caused by parity violating interaction within the nucleus, but should favourably be probed 
in atomic and molecular experiments. In atoms with stable nuclei, nuclear-spin independent 
terms caused by exchange of Z° bosons between nucleus and electrons typically dominate 
parity violating effects and often mask those effects depending on the nuclear spin that are 
significantly smaller. Thus, as of yet, only for one nucleus, namely 133 Cs, nuclear spin- 



dependent parity violating effects could be determined in atomic experiments [13j]. In linear 
open-shell molecules, the special electronic structure itself suppresses the contribution from 
nuclear spin-independent P-odd terms and offers, in principle, convenient access to nuclear- 
spin dependent V-odd contributions for a variety of nuclei, including those with an odd 
number of neutrons instead of an odd number of protons. A complementary route would be 



the detection o 
techniques \l4 - 



MST> — VV in polyatomic chiral molecules by nuclear magnetic resonance 



23]. As of yet, however, molecular parity violation has not been detected, 
which underlines the particular need for identification of promising molecular candidate 
systems by theoretical means. 

In open-shell diatomic molecules, the contribution from interactions in Eq. ([1]) to the effec- 
tive molecular spin-rotational Hamiltonian can be parametrised by the term W a k^ [XxS }-I 



24[ , a parity- violating contribution to the hyperfine coupling tensor, where A is the unit vec- 
tor pointing along the molecular axis, S eS is the effective electron spin and W a is a single 
constant characterising the P-odd electron spin-nuclear spin coupling for a given nucleus 
with nuclear spin I. In the basis of the degenerate molecular states \Q) and |— fi) (Q is the 
projection of the total electronic momentum on the molecular axis coinciding with the z 
axis) W a is approximately (see discussion below) proportional to the non-diagonal matrix 



element of the operator in Eq. (JTJ) : 

1 dh 1 

W & = (Q\—^-\-Q) xv , (3) 

where it was accounted for that A has only a no n- vanishing z-component. In contrast to W a , 
components of the hyperfine coupling tensor A can be computed also as diagonal matrix 
elements in the \Q) and |— Q) basis, which we will exploit below to estimate spin-polarisation 
effects on W a . 

To calculate W a we utilise a quasi-relativistic two-component zeroth-order regular ap- 
proximation (ZORA) approach to electroweak quantum chemistry, which proved to perform 
well in calculations of the spin-independent P-odd energy differences for chiral compounds 
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26]. Details of the ZORA approach 



and below we give only the 



when compared to four-component treatment 
for one and multielectron cases can be found elsewhere 
final expression of MSD — VV terms in the ZORA approach in the self-consistent field (SCF) 
framework of Hartree-Fock-Coulomb and Kohn-Sham-Coulomb. The derivation of these 



terms together with common consideration of the parity violation problem in open-shell 



polyatomic molecules can be found in Ref. 28] 



TABLE I: Parity violating terms in ZORA Hamiltonian. 

Term Name Expression 

Zs^' 1 ^ Scalar P-odd interaction ^^ Qa{^ 1 V-, Pa(^)} 

z hf Scalar P-odd / hyperfine P-cvcn interaction ^2-^2 Q A { e & ' ^-^(r), Pa(^)} 

2 8d Nuclear spin-dependent P-odd interaction 2^/2 ^-^i A { & ' P> ■ IaPa(t)} 

Z 8dr Nuclear spin-dependent P-odd / hyperfine P-evcn interaction ^-^=2fc^ tJ 4 {e <r " ^^{r), ^- & ■ Tjipj^(r)} 

In Table HJ Qa is the weak charge of nucleus A, Qa = Na — (1 — 4 sin 2 9 w )Za, where 
Na is the number of neutrons in nucleus A, Za the nuclear charge, sin 2 9w the Weinberg 
parameter, for which we employ the numerical value sin 2 6\v = 0.2319, and A^ is the magnetic 
vector potential from the point-like nuclear magnetic moments ft a = TvyaIa with A^(r) = 
(/xo/47r) x ~ Ra)/{\t — Ra\) 3 , 1a being the gyromagnetic ratio and /io being the 

vacuum permeability. The symbol e denotes the elementary charge (charge of a positron), 
m e the mass of the electron, h = h/ (2tt) the reduced Planck constant and {x, y} = xy + yx 
the anticommutator. The ZORA factor Co is also used, Q = 1/ ^2m e — V/jf^j, where V is the 



model potential (with additional damping 29]) proposed by van Wiillen [30|, which alleviates 



the gauge-dependence of ZORA. To calculate W a , the terms of the ZORA Hamiltonian which 
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are first order in / have to be accounted for, namely 



(0,1) + ^ = E H (W e a ■ A^l, V(r)} + 2Au,4<? • p, - c a ■ T A p A {f)}^ 



(4) 



An advantage of the ZORA approach is that one of the terms coupling the P-even hyperfine 
interaction with the P-odd nuclear spin-independent weak interaction (the first term in 
Eq. (jlj)) naturally appears after the transition from a four-component to a two-component 
framework. In our calculations we neglect this term together with accompanying response 
terms, however, as in atomic calculations it was shown to give corrections on the order of a 
few percent for heavy atoms. 



The Z-dependent scaling behaviour o 



P-odd interaction was first obtained in jl | and, for nuclear spin-dependent P-odd interaction, 
in J 

(s 1/2 |^ v |p 1/2 ) ~ Z 2 R(Z) (5) 

R(7) = 4 2 > /l-(Za)2 + l / a V"^ 1 
1 ' 3r(2Vl-(W + D 2 V MAW. r J ' W 



the matrix element of the nuclear spin-independent 



where ao is the Bohr radius, a the fine structure constant and tq = 1.2 fm a nuclear size 
parameter. The analytic form of the relativistic enhancement factor R(Z) was obtained 
from a model treatment, such that R(Z) can either be tabulated and used in approximate 
treatments or be calculated directly in atomic relativistic vs. non-relativistic calculations. As 
is shown in Fig.[T]on a double logarithmic scale, R(Z) depends heavily on Z (here the analytic 
form was used). The proposed scaling behaviour for atomic systems can subsequently be 
studied in explicit calculations for molecular systems as we will show below. 



CALCULATION DETAILS AND RESULTS 



The results of our study are summarised in Table IIHI We calculated the absolute value 
of the parameter W a for the ground £1/2 states of the alkaline earth metal monofluorides 
(Mg-Ra)F and group XII monohydrides (Zn-Cn)H. In addition, we report results for YbF 
for comparison with other approaches. In all our computations we used for the alkaline 
earth metal atoms a basis set of uncontracted Gaussians with the exponent coefficients (EC) 
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FIG. 1: Relativistic enhancement factor R(Z) as a function of the nuclear charge Z shown 

on a double logarithmic scale. 



composed as an even-tempered series. This sequence was generated according to ctj = r y[5 N ~ t , 
i = 1, . . . , N. For s,p-functions was taken equal to 2.0 for basis sets centred on the heavier 
alkaline earth metal nuclei (Sr, Ba, Ra) as well as on ytterbium and (5/2) x / 25 x 10 2 / 5 « 2.6 
for the group XII nuclei (Zn, Cd, Hg, Cn) and the lighter group II nuclei (Mg, Ca). For 
all sets of d-functions and f-functions f3 = (5/2) 1 / 25 x 10 2 / 5 ~ 2.6 was chosen. The tighter 
basis sets for s- and p-functions were employed, because the P-odd operator mixes mainly 
s- and p-waves on the heavy nucleus. For the ECs and the resulting basis set dimensions, 
see TableJTIl On the fluorine atom in all cases an uncontracted ANO basis set of triple-zeta 
quality |3l| and on hydrogen an s,p-subset of an uncontracted correlation-consistent basis 



set of quadruple-zeta quality 



3 11 ] were used with the ECs given explicitly in Table [TT1 



,R 2 



The nuclear density was modelled by a spherical Gaussian distribution p(R) = po e 25 ' 
where £ is the mean root square radius of the corresponding nucleus computed according to 
the empirical formula £ = (0.836A 1 / 3 +0.57) fin = (\.hl§%A x l z +\mi) 10" 5 a , where A is the 
given mass number of the respective isotope. Within this work we employed mass numbers 
corresponding to the standard relative atomic mass rounded to the nearest integer, namely 
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X H, 19 F, 24 Mg, 40 Ca, 65 Zn, 88 Sr, 112 Cd, 137 Ba and 201 Hg. The radium nucleus with atomic 
mass number 225 was taken and the copernicium nucleus with atomic mass number 284. As 
computed W a values do not depend too strongly on the atomic mass number (for 213 RaF, 
223 RaF and 225 RaF changes in |W a | were found on the order of a few Hz |7j]), we report only 
one value for \W a \, even though the specific isotope corresponding to the standard relative 
atomic mass may have a nuclear spin quantum number of I = 0. A modified version 



23. 25 
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32] was used for the complex generalised 



28j of the program package TURBOMOLE 
SCF (Hartree-Fock or Kohn-Sham) calculations. As spatial symmetry was not exploited, 
we also calculated the value of the projection Q of the total electron angular momentum 
on the molecular axis. In the two-component generalised Hartree-Fock (GHF) calculations 
Q = 0.5 ± 10~ 3 , in the two-component density functional theory (DFT) calculations within 
the generalised Kohn-Sham (GKS) framework Q = 0.5 ± 10~ 4 . The value of |W a | was 
calculated according to Eq. ([3]). In the ZORA calculations |0) and | —Q) are many-electron 
functions and the Lowdin formula [33] for calculations of the matrix elements between non- 
orthogonal one-determinantal (OD) functions was applied: 

(*!iwi* a > = WiWio, ( 7 ) 

k I 

where \^\) and \^2) can De either orthonormalised or non-orthonormalised OD functions, 
(ipi\w\ipk) is the matrix element of the one-electron operator w between members of the two 
sets of molecular spin-orbitals with (ipi\ipj) = Sij, (ipi\ipj) = that are occupied in the OD 
wavefunctions and ^2, respectively, and T>(k\l) is the minor of S of rank n — 1 (n is the 
number of electrons) which is obtained from the original OD wavefunctions by crossing out 
in the overlap matrix S (that has the matrix elements = (ipkl^i)) the k-th row and l-th 
column with subsequently forming the determinant of the resulting submatrix. We note in 
passing that in the direct application of the present complex GHF/GKS approach only the 
absolute value of W a is immediately accessible, whereas determination of its sign requires an 
additional symmetrisation procedure, which is for the purpose of the present study, however, 
not required. 

In calculations of the (Mg-Ra)F row two different exchange-correlation (XC) functionals 
were used in a generalised Kohn-Sham DFT framework: 1) the local density approximation 
LDA and 2) the Becke 3-parameter hybrid functional containing Becke's exchange functional 
together with Lee- Yang- Parr's (LYP) correlation functional B3LYP. This latter hybrid XC 



functional, which contains an admixture of about 20 % non-local Fock exchange, was used 



in the form employed in the Gaussian 03 program package [34| with an approximation 
(VWN3) to the correlation functional of the homogeneous electron gas. The equilibrium 
distance for all diatomic molecules was taken from experimental data, except for RaF, where 
the distance was obtained in [7J from four-component relativistic coupled cluster calculations 
in the Fock space (FS-RCC-SD) and CnH, where we used the bond length obtained in the 
two-component GHF framework. 

One can see from Table [I] two main trends when accounting for correlations by DFT: 
1) systematic increase in the value of |W a | from B3LYP to LDA XC functionals and 2) 
relative decrease in correlation contributions from 33 % for MgF to about 12 % for RaF. 
Both these dependences are consistent with previous observations and anticipations. The 
former trend was observed for parity violating energy differences between enantiomers of 
chiral molecules 



26] , the latter is also not surprising as the main contribution in this class of 
heavy-atom open-shell diatomic molecules is expected to arise from spin-polarisation effects, 
which cannot (fully) be accounted for by direct calculation of non-diagonal matrix elements 
between complex GHF wavefunctions at least for T-odd operators, for which thus results 
of essentially paired GHF quality are obtained. A discussion of the influence of symmetry 
brea king for OD wavefunctions on matrix elements of different operators can be found in 



Ref. 35j. Finally we plotted on a double logarithmic scale (Fig. |2]) instead of |W a | the values 
of \W a /R(Z)\ obtained on the GHF level against Z, as we have argued previously 
that one should correct for the relativistic enhancement factors when attempting to extract 
Z k scaling laws from quasi-relativistic and relativistic calculations. Fitting of the points in 
Fig. [2] by a linear function gives a slope equal 1.75 for (Mg-Ra)F and 2.68 for (Zn-Cn)H 
which is indeed close to the scaling factor for MST> — VV interaction matrix element. 
Our current results have been mainly confirmed by recent four-component calculations of 



W a in the series of diatomic radicals (Mg-Ra)F [36j]. The authors of [36| have also observed 
the Z k scaling (also with k close to 2) for W a /R(Z). Besides performing Dirac-Hartree- 
Fock-Coulomb and Dirac-Kohn-Sham calculations in a paired GHF and paired GKS frame- 



work, which can not account for core-polarisation effects, the authors of Ref. [36] employed 
some approximate atom based schemes to roughly estimate part of electron correlation ef- 
fects via scaling factors. Our treatment, however, is based on the complex GHF/GKS 
framework and thus allows to capture part of the electron correlation effects directly within 



the molecular calculations (see also below), whereas some contributions are not included 
due to calculation via off-diagonal matrix elements between time-reversed wavef unctions. 
Our direct DFT-based estimates for \W a \ in MgF, CaF and SrF can reasonably well be 
reproduced by the indirect procedure employed in Ref . [3^ . For BaF and RaF we find, how- 
ever, only a modest electron correlation contribution on the DFT level of theory whereas in 
Ref. jsfj significant changes are reported for RaF. Even without the subsequent attempts 
to account for further electron correlation effects, in Ref. [36] DHF and DKS values for |W a | 
in RaF differ by more than 15 % (and by about 30 % for LDA XC functional) and in BaF 
only by 2 %. As it was ment.oned the latter result (for BaF) is actually iu agreement with 
the earlier calculations in Ref. [37J, in which the authors found that the main contribution 
comes from the spin-polarisation effects, though further accounting for electron correlation 
gives minor contribution. We note in passing that in Ref. [36| the data are mixed up for 
the calculations without accounting for spin-polarisation (SCF in notations of the authors of 
Ref. 0, W a = 111 Hz), with the accounting for spin-polarisation (SCF-EO, W a = 181 Hz) 
and electron correlation + spin-polarisation (RASSCF-EO, W a = 175 Hz). In Ref. 36 ] 
Faegri's energy-optimised basis sets were employed, which required augmenting with addi- 
tional functions to be used in calculations of properties that depend on the behaviour of 
the electronic wavefunction near the nucleus. To check the influence of the basis set choice, 
we performed calculations with two additional basis sets (see below for basis set specifica- 
tion) for Ra together with an uncontracted aug-cc-pVTZ basis set Ils6p3d2f on fluorine 
nucleus. The first Ra basis set (Basis S in Table IIIII) was Faegri's uncontracted basis set 



25s21pl4d9f recommended for relativistic calculations 38j and another one (Basis L) was a 



large even-tempered basis set 36s33p22dl5f generated according to recommendation of the 



article 39|. The result of the calculations with these basis sets clearly shows that with the 
extension of the basis set from Basis S to Basis L the difference between GHF and GKS 
results of essentially paired generalised SCF quality decreases from 38 % to 20 % for LDA 
XC functional, getting close to the values reported by us (~ 12 %) for the basis sets we 
used herein and in Ref. [3]. This prov ides some indication that the pronounced electron 



correlation effects reported in Ref. 36( for the DFT framework might primarily be caused 
by the special choice of basis set therein. 

To estimate spin-polarisation contributions within the GHF approach we use scaling re- 



lations from the semiempirical molecular model by Kozlov 



40j, which is known to reproduce 



ab initio parameters of the V, T-odd spin-rotational Hamiltonian for ground states of BaF, 
YbF and some other molecules with an accuracy of 10-15%. For this model some simple (ap- 
proximate) relations can be established between the parameters of the electronic structure, 
required for calculations of W a of the linear diatomic molecules employed in the current work 
in their E1/2 ground states, and the hyperfme coupling tensor terms A iso (isotropic) and A d 
(dipole). For our purpose and the current set of molecules, however, more important is that 
the relation between W a obtained for different approximations (e.g. complex and paired 
generalised Hartree-Fock wavef unctions, cGHF and pGHF, respectively) is approximately 
equal to the ratio between the square root of the products of A iso and A d (see equations 



(33), (34) and (10) in Ref. [40|; assuming that signs of A iso , A d are identical): 

(8) 



^ cGHF r(Aso-A d ) cGHF ^ 1/2 



Wi GilF l(Ai SO • Ai) pGHF J ' 
Thus, by calculating the hyperfme tensors with accounting for spin-polarisation (in our case 
as diagonal matrix elements within the complex GHF scheme) and without accounting for 
it (as non-diagonal matrix elements, leading to results of essentially paired GHF quality) 
we can restore spin-polarisation contributions, which are expected to be most important for 
RaF. The results of this scaling are presented in Table HVl One can observe that for molecules 
with a valence electronic structure similar to RaF the relative deviation of our scaling for W a 
parameters is better than 10 % when judged from the corresponding RASSCF/EO results 



or about 10 % in comparison with the semiempirical estimates for HgH by Kozlov 40] . This 
finding is particularly encouraging for the identification of promising molecular candidates, 
although we expect the accuracy in general to be somewhat lower than implied by the present 
results. Thus, for reliable estimates of electron correlation and spin-polarisation effects on 
the value of W a (and other properties depending on the behaviour of the wavefunction 
near the nucleus) in RaF one has to employ high-order correlation calculations, for instance 



similar to those in Ref. 1411. It is also interesting to note that in the group XII monohydrides 
series the spin-polarisation contribution should suppress rather than enhance the MST> — 
VV interaction, at least for the two heavy representatives reported in Table IIVI 

Finally, we emphasise that although the treatment of relativistic effects in the four- and 
two-component framework is different, deviations between results for MST> — VV operators 
are not expected to be significantly larger than 3 % for the heavier nuclei (row 4-7) when 



judged on the basis of earlier calculations 



23 
10 



26|, provided appropriate basis sets are used 



(see also discussion in Ref. [3j]). 




FIG. 2: Scaling of \WJ R(Z)\ GHF values (in Hz) with Z for the (Mg-Ra)F and (Zn-Cn)H 
series (orange line, blue dots and red line, green dots, respectively) on a double logarithmic 
scale. The slope of the lines is 1.75 for (Mg-Ra)F and 2.68 for (Zn-Cn)H, which implies a 
R(Z)Z k scaling law for W a with k = 1.75 and k = 2.68, respectively. 
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TABLE II: Basis sets parameters for ZORA HF/DFT calculations. Even-tempered basis 
sets of uncontracted Gaussians are given in the form N^; I; (EC max ; EC m ; n ), where iV"b as is 
the number of Gaussians, / is s, p, d or f and represents the angular momentum quantum 
numbers 0, 1, 2 or 3. EC max and EC min are the largest and smallest exponent coefficients, 

respectively. 





Mg 




Ca 


Sr,Ba,Yb 




Ra 


27; s; 


(500000000; 0.00769) 


27; s; 


(500000000; 0.00769) 37; s; 


(2000000000; 0.0291) 


39; s; 


(2000000000; 0.00728) 


25; p; 


(191890027; 0.02000) 


25; p; 


(191890027; 0.0200) 34; p; 


( 500000000; 0.0582) 


34; p; 


( 500000000; 0.0582) 


4; d; 


(3.750; 0.21336) 


13; d; 


(13300.758; 0.135789) 14; d; 


(13300.758; 0.0521) 


14; d; 


(13300.758; 0.0521) 








9; f; 


(751.8368350; 0.3546) 


10; f; 


(751.8368350; 0.13638) 



27; s 
25; p 
14; d 
8; f; 



Zn,Cd,Hg,Cn 
(500000000; 0.0077) 
(191890027; 0.0200) 
(13300.758; 0.0521) 
(751.8368350; 0.9219352) 



103109.46 
15281.007 
3441.5392 
967.09483 
314.03534 
113.44230 
44.644727 
18.942874 
8.5327430 
3.9194010 
1.5681570 
0.6232900 
0.2408610 
0.0843010 



F ANO basis 
P 

245.33029 
56.919005 
17.604568 
6.2749950 
2.4470300 
0.9950600 
0.4039730 
0.1548100 
0.0541840 



d 

5.000000 
1.750000 
0.612500 
0.214375 



82.640 

12.410 
2.8240 
0.7977 
0.2581 
0.08989 



2.2920000 
0.8380000 
0.2920000 
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TABLE III: Calculated "P-odd parameter \W a \ (in Hz) for open-shell diatomic molecules 
together with the charge number Z of the heavy nucleus and the equilibrium distance R e 

employed. 



t^al/Hz 







Z R e /ao 


GHF 


GKS/B3LYP 


GKS/LDA 


MgF 




12 


3.30 


3.9 


4.9 


5.2 


CaF 




20 


3.71 


8.0 


9.2 


9.5 


SrF 




38 


3.92 


3.9 x 10 1 


4.6 x 10 1 


4.8 x 10 1 


BaF 




56 


4.07 


1.11 x 10 2a 


1.19 x 10 2 


1.25 x 10 2 


RaF 




88 


4.24 


1.30 x 10 3b 


1.42 x 10 3 


1.47 x 10 3 


RaF (Basis 


S) 


88 


4.24 


1.07 x 10 3 




1.48 x 10 3 


RaF (Basis 


L) 


88 


4.24 


1.24 x 10 3 




1.50 x 10 3 


ZnH 




30 


3.01 


4.7 x 10 1 






CdH 




48 


3.36 


2.23 x 10 2 






HgH 




80 


3.33 


3.30 x 10 3c 






CnH 




112 


3.10 


4.88 x 10 4 







- 



a) In Ref. 7] 111 Hz were obtained with a slightly different basis set. 

b) In Ref. M] 1.3 kHz were reported for a slightly different basis set. 

c) In Ref. [7J 3.3 kHz were reported for a slightly smaller basis set. 
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TABLE IV: Calculated and scaled ab initio values for the parameter |W a | in BaF, YbF 

1/2 



and RaF together with scaling factor / 



(A iso -A d )^ HF 



Additional ab initio results and 



calculation methods are taken from the corresponding references. A bond length of 3.80 

was used in the calculation for YbF. 



|wg/Hz 



Method 



Ref. 



GHF 



Scaled 



ab initio 



f 



BaF 

YbF 
RaF 
HgH 
CnH 



1.11 x 10 2 

4.65 x 10 2 
1.30 x 10 3 
3.30 x 10 3 
4.88 x 10 4 



1.9 x 10 2 

6.1 x 10 2 
2.1 x 10 3 

2.0 x 10 3 

3.1 x 10 4 



1.81 x 10 2 1.68 SCF/EO 
1.75 x 10 2 RASSCF/EO 
6.34 x 10 2 1.31 RASSCF/EO 
1.65 

a 0.62 
0.63 



[22 
W 

[42] 



a) Semi-empirical estimate of Ref. 



4Cj based on spectroscopic parameters of 199 HgH and 201 HgH 



W a = 1800 Hz and W a = 1940 Hz, respectively. 
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CONCLUSIONS 



We have reported herein a numerical study on nuclear charge dependent scaling of molec- 
ular properties in open-shell diatomic molecules. After we have accounted for a relativistic 
enhancement factor R(Z), which grows non-polynominally with the nuclear charge Z, we 
obtain an approximate Z 2 scaling behaviour for the nuclear spin-dependent parity violating 
parameter W a computed at the respective equilibrium structures. This term contributes to 
the effective spin-rotational Hamiltonian used for high-resolution studies which aim for the 
first detection of molecular parity violation. The present confirmation of a simple scaling 
law is excellent news as it allows for quick estimates of parity violating effects in a whole 
series of diatomic molecules. Within the complex generalised Kohn-Sham framework em- 
ployed in this work, part of electron correlation effects on this molecular property can be 
accounted for, although some contributions are still missing. Spin-polarisation effects can 
approximately be included within a simple, but powerful molecular scaling scheme utilised 
previously for semi-empirical estimates. For high accuracy calculations, more sophisticated 
molecular electron correlation approaches are clearly needed, but given the present exper- 
imental status, the current approximate approaches allow to identify promising molecular 
candidates such as RaF, which was proposed in Ref. jz|. 
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